Cosmic magnetization: from spontaneously emitted aperiodic turbulent to ordered 

equipartition fields 
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It is shown that an unmagnetized nonrelativistic thermal electron-proton plasma spontaneously 
emits aperiodic turbulent magnetic field fluctuations of strength \SB\ = 9j3 e g 1 ^ 3 Wl^ 2 G, where 
/3 e is the normalized thermal electron temperature, W e the thermal plasma energy density and g 
the plasma parameter. Aperiodic modes fluctuate only in space, but are not propagating. For 
the unmagnetized intergalactic medium, immediately after the reionization onset, the field strength 
from this mechanism is about 4.7 ■ 10~ 16 G, too weak to affect the dynamics of the plasma. The 
shear and/or compression of the intergalactic medium exerted by the first supernova explosions 
amplify these seed fields and make them anisotropic, until the magnetic restoring forces affect the 
gas dynamics at ordered plasma betas near unity. 

PACS numbers: 07.55.Db, 52.25.Gj, 52.35.Ra, 95.30.Qd, 98.62. Ra 

Keywords: generation of magnetic fields - fluctuation phenomena - plasma turbulence - intergalactic matter 
- astrophysical plasmas 



The interstellar medium (ISM) is filled with (1) a dilute 
mixture of charged particles, atoms, molecules and dust 
grains, referred to as interstellar gas and dust, (2) par- 
tially turbulent magnetic fields, (3) dilute photon radia- 
tion fields from stars, dust and the universal microwave 
background radiation, and (4) cosmic ray particles with 
relativistic energies. It is known for a long timei~— , even 
before the discovery of the universal cosmic microwave 
background radiation, that these ISM components have 
comparable energy densities and pressures, each of the 
order of 10 -12 erg cm~ 3 , commonly referred to as the 
equipartition condition in the ISM. Since today, this truly 
remarkable equipartition in the ISM has not been under- 
stood nor explained theoretically. One refers to pressure 
partition, if the ratio of any two individual pressures is 
a constant, and to pressure equipartition, if this ratio is 
near unity. 

In other astrophysical objects equipartition conditions, 
and the closely related minimum-energy assumption, for 
the total magnetic field energy density and the kinetic 
energy density of plasma particles are also often invoked 
for convenience^ in order to analyze cosmic synchrotron 
intensities. The minimum-energy assumption was first 
proposed by BurbidgeS. and applied to the optical syn- 
chrotron emission of the jet in M87. Duric£ argued that 
any major deviation from equipartition would be in con- 
flict with radio observations of spiral galaxies. Observa- 
tionally, for a variety of nonthermal sources the equipar- 
tition concept is supported by magnetic field estimates as 
e.g. in the Coma cluster of galaxies and radio-quiet active 
galactic nuclei^. Also the solar wind plasma exhibits near 
equipartition conditions: ten years of Wind/SWE satel- 
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lite data£ near 1 AU showed that the proton and electron 
temperature anisotropies A — T±/Tu are bounded by the 
mirror and firehose instabilities at large values of the par- 
allel plasma beta /3y = StrnksT^/ B 2 > 1 and by parallel 
propagating Alfven waves 9 at small values of /3|| < 1, 
resulting in near magnetic field equipartition. 

Because of their comparably low gas densities, all 
cosmic fully and partially ionized non-stellar plasmas 
are collision-poor, as indicated by the very small val- 
ues of the plasma parameter g — v ee /iu p ^ e < 10 -10 , 
given by the ratio of the electron-electron Coulomb col- 
lision frequency v ee to the electron plasma frequency 
u>p, e , characterizing interactions with electromagnetic 
fields, so that fully kinetic plasma descriptions are nec- 
essary. Because of the large sizes of astrophysical sys- 
tems compared to the plasma Debye length, the fluctu- 
ations are descibed by real wave vectors (k) and com- 
plex frequencies u)(k) — r(k) + i'j(k), implying for the 
space- and time-dependence of e.g. magnetic fluctua- 
tions the superposition of S B{x, t) cx exp[i(fc-a?— rt)+ji\. 
One distinguishes between collective modes with a fixed 
frequency-wavenumber dispersion relation, also referred 
to as normal modes, and non-collective (no frequency- 
wavenumber relation) modes in the system. Regarding 
frequency, basically two fundamental types of fluctua- 
tions occur: (1) weakly amplified/damped solutions (e.g. 
Alfven waves, electromagnetic waves) with |7| <C wjj, 
and (2) weakly propagating solutions (e.g. firehose and 
mirror fluctuations) with ojr <C 7, including aperiodic 
solutions with ojr = (e.g. Weibel fluctuations^). Ape- 
riodic modes fluctuate only in space, do not propagate as 
r = 0, but permanently grow or decrease in time depend- 
ing on the sign of 7. Past research-^— has concentrated 
predominantly on the fluctuations from collective weakly 
amplified modes in the plasma. 

All plasmas, including unmagnetized and those in ther- 
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mal equilibrium, have fluctuations so that their state 
variables such as density, pressure and electromagnetic 
fields fluctuate in position and time. Unlike for weakly 
amplified/damped modes, however, for aperiodic fluc- 
tuations the expected fluctuation level has never been 
calculated quantitatively. Only recently general expres- 
sions for the electromagnetic fluctuation spectra (elec- 
tric and magnetic field, charge and current densities) 
from uncorrelated plasma particles in unmagnetized plas- 
mas for arbitrary frequencies have been derived^ using 
the system of the Klimontovich and Maxwell equations, 
which are appropriate for fluctuations wavelengths longer 
than the mean distance between plasma particles, i.e. 

k < fcmax = 27m;!/ 3 . These general expressions are co- 
variantly correct within the theory of special relativity, 
and hold for arbitrary momentum dependences of the 
plasma particle distribution functions and for collective 
and non-collective fluctuations. The electric^ and mag- 
netic field fluctuations in unmagnetized plasmas with the 
plasma frequency uj 2 a = (4Tre 2 n a /m a ) 1 / 2 
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are given in terms of the parallel and perpendicular form 
factors 
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and the general longitudinal and transverse dispersion 
functions Ai j r(A;,w) ivolving the respective parallel and 
perpendicular dielectric tensor elements. Eqs. ([J) - ((2]) 
are the generalizations of the standard expression found 
in the literature in which the weak, amplification limit of 
7 — >• + is taken at the outset to approximate the factor 
— i(j + i(k-v — r)) _1 by lim 7 ^ +(~ i)("f + t(k-v — r)) _1 — > 
it 5(k ■ v — r). 

We now consider the unmagnetized intergalactic 
medium (IGM) immediately after the reionization onset, 
assuming that any earlier cosmological magnetization has 
vanished during the long recombination era with a fully 
neutral IGM. Modeling the photoionization by the first 
forming star o i 18 indicates IGM temperatures of about 

n e = 10 

thermal IGM proton-electron plasma we follow recent 
workiii to calculate from Eqs. ([T])- Eqs. ([2]) the spon- 
taneously emitted magnetic field fluctuation spectrum of 
aperiodic (r = 0) fluctuations 



T p = T = 10 T4K and ionized gas densities of 
~ 7 ri7 cm~ 3 at redshift z = 4. For this isotropic 
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with the thermal velocity u a = y // 2kBT a /m a and D(x) = 
e x erfc (\x\) denoting the complimentary error function. 
The related collective Weibel mode^ has a positive 
growth rate in anisotropic plasma distribution functions, 
but is not excited in isotropic plasma distributions. 

Integrating over all values of 7 and k provides the en- 
ergy density of spontaneously emitted fully random mag- 
netic fluctuations 
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where 



F{x,n) = D(x) +[i- 1 D{xfi), (6) 

the mass ratio /i 2 = m p /m e = 1836 and /3 e — u e /c — 

(2) 1.84- 10- 3 T 4 1/2 . Likewise, one finds for aperiodic fluctu- 
ations no charge density and parallel electric field fluctu- 



ations < Sp 2 >k=< SE 2 > fc = 0, and 
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The symmetric function D(—x) = D(x) attains its max- 
imum value -D(O) = 1 at x = and is monotonically 
decreasing to as x — > 00. Its rational approximation 
better than 2.5T0 -5 is given byi2 D(x) ~ ait — a 2 i+a.3i 3 
with t = /(l + px), p = 0.47047, ai = 0.3480242, 
a 2 = 0.0958798 and a 3 = 0.7478556. Given the small- 
ness of a 2 we use as lower and upper limits Dl{x) < 
D(x) < Du{x) with Dl{x) = (ai — a 2 )t + a 3 t 3 and 
Du(x) = ait + (a 3 - a 2 )t 3 , i.e. 



D L {x) ~ (ao + a 3 t 2 )t, D v {x) ~ ( ai + a 4 t 2 )t (8) 



(4tt 5 / 2 ) < SB 2 > k , 



with ao = ai — a 2 = 0.2521444 and 04 — a 3 — a 2 = 
0.6519758. Fig. 1 indicates that the relative deviation of 
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FIG. 1: Relative variations of the lower (Dl(x)/D(x)) (dotted 
curve) and upper (Du(x)/D(x)) (dashed curve) limits as a 
function of x. 



the upper and lower limits are smaller than 30 percent 
at all values of x. 

Because of the large value of fi = 43 we can neglect 
the proton contribution in Eq. (|5|) so that F(x, fi) ~ 
D(x), implying in terms of the normalized wave vector 
k = kc/u) Pt e that 
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Eq. Q then becomes 
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The integral (fT0|) is well approximated by Jl(/3,k) < 
J(f3, k) < Ju{P, k) with 
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After straigtforward but tedious algebra we derive 
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with Y^ u( K ) — iT^^ax o/\pn 2 ]. The asymptotic expan- 
sions for small and large values of Yl,c/ correspond to 
large and small values of the normalized wavenumber, 
respectively, providing 
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The resulting expressions for < 5B > k ' from Eq. (J5J) 
to leading order increase oc « at small normalized wave- 
length k < (aifiit 1 / 2 /p) 1 / 2 and approach constants at 
large k, disagreeing with earlier result a 20 ' 21 which were 
based on fluctuation spectrum formula valid for weakly 
damped/amplified modes. The constants at large values 
of k provide the dominating contribution to the remain- 
ing K-integral in Eq. ([IT]). We find < SB 2 > v = 1.38 < 
SB 2 >l with 
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with the thermal energy density W e = 3n e kBT/2 = 
2.1 ■ 10~ 19 njTx erg cm -3 and the plasma beta g = 
2.3 ■ 10 _13 n^ 2 T 4 3 ^ 2 . This magnetic field energy den- 
sity corresponds to a minimum total fluctuating magnetic 
field strength of 



\SB\ 



< SB >, 



,1/2 
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This fluctuating magnetic field strength is still small 
enough to allow rectilineal propagation of ultrahigh en- 
ergy cosmic ray protons in the IGM, The associated tur- 
bulent plasma beta is larger than 



87m e fcsT 1 

< SB 2 > v = 6.9^2/3 



= 1.1 • 10 
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Now we finally address how ordered magnetic field struc- 
tures emerge from these randomly distributed magnetic 
fields. As we demonstrated the unmagnetized, isotropic, 
thermal and steady IGM plasma by spontaneous emission 
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generates steady tangled fields, isotropically distributed 
in direction, on small spatial scales < 10 10 n 7 ^ 2 cm (cor- 
responding to k > 1). Because of its ultrahigh turbulent 
plasma beta value ([T7|) . these seed helds are too weak 
to affect the dynamics of the IGM plasma, but are tied 
passively to the highly conducting IGM plasma as frozen- 
in magnetic fluxes. Earlier analytical considerations and 
numerical simulations^— showed that any shear and/or 
compression of the IGM medium enormously amplify 
these seed magnetic fields and make them anisotropic. 
Considering a cube containing an initially isotropic mag- 
netic field, which is compressed to a factor 77 <C 1 times its 
original length along one axis (z), these authors showed 
that the perpendicular magnetic field components are en- 
hanced by the factor r; -1 . Depending on the specific 
exerted compression and/or shear, even one-dimensional 
ordered magnetic field structures can be generated out 
of the original isotropically tangled field configuration 26 . 
Hydrodynamical compression or shearing of the IGM 
medium arises from the shock waves of the supernova 
explosions of the first stars at the end of their lifetime, 
or from supersonic stellar and galactic winds. The IGM 
seed magnetic field upstream of these shocks is random 
in direction, and by solving the hydrodynamical shock 
structure equations for oblique and conical shocks it has 
been demonstrate d 27 ! 28 , that the shock compression en- 
hances the downstream magnetic field component paral- 



lel to the shock, but leaving the magnetic field component 
normal to the shock unaltered. Consequently, a more or- 
dered downstream magnetic field structure results from 
the randomly oriented upstream field. Such stretching 
and ordering of initially turbulent magnetic fields is also 
seen in the numerical hydrodynamical simulations of su- 
personic jets in radio galaxies and quasars^. 

This passive hydrodynamical amplification and 
stretching of magnetic fields continues until the mag- 
netic restoring forces affect the gas dynamics, i.e. at 
ordered plasma betas near unity. During these stretch- 
ings and amplifications also MHD dynamo processes 
need to be considered 2 ^. The stretching and ordering 
of magnetic fields also affects the motions of the IGM 
protons and electrons, so that our original expressions 
for the spontaneously emitted fluctuation spectra ([1]) in 
unmagnetized plasmas then no longer apply. 

In principle, our suggested mechanism of sponta- 
neously emitted aperiodic turbulent magnetic fields 
should also operate during earlier cosmological epochs 
before recombination. However, the then relativistic 
temperatures of the unmagnetized plasma require the 
use of the relativistic fluctuation spectra which also are 
available^. 
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